LONG TIME EXISTENCE OF SMOOTH SOLUTIONS FOR 
THE RAPIDLY ROTATING SHALLOW- WATER AND EULER EQUATIONS 



BIN CHENG AND EITAN TADMOR 

Abstract. We study the stabilizing effect of rotational forcing in the nonlinear setting of two- 
dimensional shallow-water and more general models of compressible Euler equations. In [T^ we have 
shown that the pressureless version of these equations admit global smooth solution for a large set of 
sub-critical initial configurations. In the present work we prove that when rotational force dominates 
the pressure, it prolongs the life-span of smooth solutions for t ^ ln(5~^); here 5 ^ 1 is the ratio of 
the pressure gradient measured by the inverse squared Fronde number, relative to the dominant rota- 
tional forces measured by the inverse Rossby number. Our study reveals a "nearby" periodic-in-time 
approximate solution in the small 5 regime, upon which hinges the long time existence of the exact 
smooth solution. These results are in agreement with the close-to periodic dynamics observed in the 
"near inertial oscillation" (NIO) regime which follows oceanic storms. Indeed, our results indicate the 
existence of smooth, "approximate periodic" solution for a time period of days, which is the relevant 
time period found in NIO obesrvations. 
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1. Introduction and statement of main results 

We are concerned here with two-dimensional systems of nonlinear Eulerian equations driven by 
pressure and rotational forces. It is well-known that in the absence of rotation, these equations 
experience a finite-time breakdown: for generic smooth initial conditions, the corresponding solutions 
will lose C^-smoothness due to shock formation. The presence of rotational forces, however, has 
a stabilizing effect. In particular, the pressureless version of these equations admit global smooth 
solutions for a large set of so-called sub-critical initial configurations, [TT]. It is therefore a natural 
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extension to investigate the balance between the regularizing effects of rotation vs. the tendency of 
pressure to enforce finite-time breakdown (we mention in passing the recent work [21j on a similar 
regularizing balance of different competing forces in the ID Euler-Poisson equations). In this paper 
we prove the long-time existence of rapidly rotating flows characterized by "near-by" periodic flows. 
Thus, rotation prolongs the life-span of smooth solutions over increasingly long time periods, which 
grow longer as the rotation forces become more dominant over pressure. 

Our model problem is the Rotational Shallow Water (RSW) equations. This system of equations 
models large scale geophysical motions in a thin layer of fluid under the influence of the Coriolis 
rotational forcing, (e.g. [El §3.3], [13 §2.1]), 

(1.1a) dth + V-{hu) = 0, 

(1.1b) 9iU + u • Vu + gV/i - /u-^ = 0. 

It governs the unknown velocity field u := {u'^^\t,x,y)^u^'^\t,x^y)^ an height h := h{t,x,y), where 
g and / stand for the gravitational constant and the Coriolis frequency. Recall that equation (jl.lap 
observes the conservation of mass and equations (jl.lbp describe balance of momentum by the pressure 
gradient, ^V/i, and rotational forcing, /u-*- := f{vP'\ — u^"*^-'). 

For convenience, we rewrite the system (jl.ip in terms of rescaled, nondimensional variables. To 
this end, we introduce the characteristic scales, H for total height /i, D for height fluctuation h — H, 
U for velocity u, L for spatial length and correspondingly L/U for time, and we make the change of 
variables 

u = u' , x'L, y'Lj U, h = H + h'(^^, x'L, y'L^ D. 
Discarding all the primes, we arrive at a nondimensional system, 

dth + u ■ Vh + + V ■ u = 0, 

dtu + u-Vu + ^Vh- = 0. 

We are concerned here with the regime where the pressure gradient and compressibility are of the 
same order, ^—r ~ — . Thus we arrive at the (symmetrizable) RSW system, 

(1.3a) dth + u-Vh+ (^^ + h^V -u = 0, 

1 1 

(1.3b) aiU + u-Vu + -V/i--Ju = 0,. 

a T 

Here a and r, given by 

{1.3c) " " 



are respectively, the Froude number measuring the inverse pressure forcing and the Rossby number 
measuring the inverse rotational forcing. We use J to denote the 2x2 rotation matrix J := 



1 

-1 

To trace their long-time behavior, we approximate (jl.3ap . (jl.3bp with the successive iterations, 

(1.4a) dthj + Uj_i • Vhj + (^^ + ^ " "i-i = 0' i = 2, 3, . . . 

(1.4b) dtUj + Uj -Vuj + -Vhj - -JvLj = 0, j = l,2, 
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subject to initial conditions, hj{0,-) = /lo(') and Uj(0,-) = uo(-). Observe that, given j, (jl.4p are 
only weakly coupled through the dependence of Uj on hj, so that we only need to specify the initial 
height hi- Moreover, for a ^ t, the momentum equations ()1.3bp are "approximately decoupled" 
from the mass equation ()1.3ap since rotational forcing is substantially dominant over pressure forcing. 
Therefore, a first approximation of constant height function will enforce this decoupling, serving as 
the starting point of the above iterative scheme. 



(1.4c) hi = constant. 

This, in turn, leads to the first approximate velocity field, ui, satisfying the pressureless equations, 

(1.5) atUi + ui-Vui-ijui = 0, ui(0,-) =uo(-). 

r 

Liu and Tadmor [17] have shown that there is a "large set" of so-called sub-critical initial configurations 
uq, for which the pressureless equations (jl.Sp admit global smooth solutions. Moreover, the pressureless 
velocity ui(t, •) is in fact 27rr-periodic in time. The regularity of ui is discussed in Section [5J 

Having the pressureless solution, (hi = constant, ui) as a first approximation for the RSW solution 
(/i, u), in Section [3] we introduce an improved approximation of the RSW equations, (/i2>U2), which 
solves an "adapted" version of the second iteration (j = 2) of ()1.4p . This improved approximation 
satisfies a specific linearization of the RSW equations around the pressureless velocity ui, with only 
a one-way coupling between the momentum and the mass equations. Building on the regularity and 
periodicity of the pressureless velocity ui, we show that the solution of this linearized system subject 
to sub-critical initial data (/io,uo), is globally smooth; in fact, both h2{t,-) and U2(t, •) retain 27rT- 
periodicity in time. 

Next, we turn to estimate the deviation between the solution of the linearized RSW system, (/12, U2), 
and the solution of the full RSW system, (/i,u). To this end, we introduce a new non-dimensional 
parameter 

5:=^ 



a2 fLU' 

measuring the relative strength of rotation vs. the pressure forcing, and we assume that rotation is 
the dominant forcing in the sense that 5 <C 1. Using the standard energy method we show in Theorem 
14. II and its corollary that, starting with sub-critical initial data, the RSW solution [h{t, •),u(t, •)) 
remains sufficiently close to (/i2(i, •)> U2(t, •)) ™ the sense that, 

\\h{t, •) - h2{t, OIIh^-S + l|u(t, •) - U2{t, ■)\\Hm-S < _ ^Cot J)2 ' 

where constant Co = Co(m, |Vuo|oo, l^oloo) • ||uo, /lollm- In particular, we conclude that for a large set 
of sub-critical initial data, the RSW equations (jl.Sp admit smooth, "approximate periodic" solutions 
for long time, t < ts '■= ln{S~^), in the rotationally dominant regime 5 <C 1. 

We comment that our formal notion of "approximate periodicity" emphasizes the existence of a 
periodic approximation (/i2,U2) nearby the actual flow (/i, u), with an up-to 0{5) ^ 1-error for 
sufficiently long time. Therefore, strong rotation stabilizes the flow by imposing on it approximate 
periodicity, which in turn postpones finite time breakdown of classical solutions for a long time. A 
convincing example is provided by the so called "near-inertial oscillation" (NIO) regime, which is 
observed during the days that follow oceanic storms, e.g. [22]. These NIOs are triggered when 
storms pass by (large U's) and only a thin layer of the oceans is reactive (small aspect ratio H/L), 
corresponding to 6 = jjjj <C 1. Specifically, with Rossby number r ~ 0(0.1) and Froude number 
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a ~ C(l) we find 6 ~ 0.1, which yield the existence of smooth, "approximate periodic" solution 
for i ~ 2 days. We note that the clockwise rotation of cyclonic storms on the Northern Hemisphere 
produce negative vorticity which is a preferred scenario of the sub-critical condition (|2.1bp . Our results 
are consistent with the observations regarding the stability and approximate periodicity of the NIO 
regime. 

Next, we generalize our result to Euler systems describing the isentropic gasdynamics, in Section 
14.21 ^-iid ideal gasdynamics, in Section [4.31 We regard these two systems as successive generalizations 
of the RSW system using the following formalism, 

(1.6a) 5tp + V-(pu) = 0, 

(1.6b) 5tu + u- Vu + /9~^Vp(/j,5) = /Ju, 

(1.6c) dtS + u-VS = 0. 

Here, the physical variables p, 5 are respectively the density and entropy. We use p(p, 5) for the 
gas-specific pressure law relating pressure to density and entropy. For the ideal gasdynamics, the 
pressure law is given as p := Ap^e^ where A, 7 are two gas-specific physical constants. The isentropic 
gas equations correspond to constant 5", for which the entropy equation (|1.6cp becomes redundant. 
Setting ^ = g', 7 = 2 yields the RSW equations with p playing the same role as h. 

The general Euler system <\1.Q\) can be symmetrized by introducing a "normalized" pressure function, 

p:=^p^{p,S), 
7-1 

and by replacing the density equation (jl.6ap with a pressure equation, 

^ - 1 

(1.6d) + u • Vp + ^^pV • u = 0. 

We then nondimensionalize the above system ()1.6bp . (jl.6cp and ()1.6dp into 

-i- 

2 \a 



7 - 1 / 1 

dtp + u • Vp H — — h p 1 V • u = 0, 



Stu + u- Vu+ e'^^Vp = iju, 



dtS + M-VS = 0. 

The same methodology introduced for the RSW equations still applies to the more general Euler 
system, independent of the pressure law. In particular, our first approximation, the pressureless 
system, remains the same as in (jl.Sp since it ignores any effect of pressure. We then obtain the second 
approximation (p2)U2,52) (or (p2)U2) in the isentropic case) from a specific linearization around the 
pressureless velocity ui. Thanks to the fact that /i, p and S share a similar role as passive scalars 
transported by u, the same regularity and periodicity argument can be employed for (p2,U2,52) 
in these general cases as for (/i2,U2) in the RSW case. The energy estimate, however, needs careful 
modification for the ideal gas equations due to additional nonlinearity. Finally, we conclude in Theorem 
4.2 and 4.3 that, in the rotationally dominant regime 5 <C 1, the exact solution stays "close" to the 
globally smooth, 27rr-periodic approximate solution (p2, U2, S2) for long time in the sense that, starting 
with H"'' sub-critical data, the following estimate holds true for time t ^ ln{6^^), 

\\pit, •) - P2it, OIU-S + ||u(t, •) - U2(t, Ollm-S + \\Sit, •) - 52(t, OIU-S < Y^^C^H' 

Our results confirm the stabilization effect of rotation in the nonlinear setting, when it interacts with 
the slow components of the system, which otherwise tend to destabilize of the dynamics. The study 
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of such interaction is essential to the understanding of rotating dynamics, primarily to geophysical 
flows. We can mention only few works from the vast literature available on this topic, and we refer 
the reader to the recent book of Chemin et. al., [6] and the references therein, for a state-of-the art of 
the mathematical theory for rapidly rotating flows. Embid and Majda [Tl|8] studied the singular limit 
of RSW equations under the two regimes r^^ ~ a^^ oo and 0(1), a ^ ^ 00. Extensions to 

more general skew-symmetric perturbations can be found in the work of Gallagher, e.g. [9]. The series 
of works of Babin, Mahalov and Nicolaenko, consult [H [H O U [5] and references therein, establish 
long term stability effects of the rapidly rotating 3D Euler, Navier-Stokes and primitive equations. 
Finally, we mention the work of Zeitlin, Reznik and Ben Jelloul in [23\ [23] which categorizes several 
relevant scaling regimes and correspondingly, derives formal asymptotics in the nonlinear setting. 

We comment here that the approach pursued in the above literature relies on identifying the limiting 
system as r ^ 0, which filters out fast scales. The full system is then approximated to a first order, 
by this slowly evolving limiting system. A rigorous mathematical foundation along these lines was 
developed by Schochet [19j, which can be traced back to the earlier works of Klainerman and Majda 
|13j and Kreiss [13] (see also [ID]). The key point was the separation of (linear) fast oscillations from 
the slow scales. The novelty of our approach, inspired by the critical threshold phenomena [16], is to 
adopt the rapidly oscillating and fully nonlinear pressureless system as a first approximation and then 
consider the full system as a perturbation of this fast scale. This enables us to preserve both slow and 
fast dynamics, and especially, the rotation-induced time periodicity. 

2. First approximation- the pressureless system 
We consider the pressureless system 

(2.1a) dtui + ui • Vui — -Jui = 0, 

r 

subject to initial condition ui(0, •) = uo(-). We begin by recalling the main theorem in |17j regarding 
the global regularity of the pressureless equations (j2.1ap . 

Theorem 2.1. Consider the pressureless equations i2.1a\) subject to C^-initial data ui(0, •) = uq{-). 
Then, the solution ui(t, •) stays for all time if and only if the initial data satisfy the critical 
threshold condition, 

(2.1b) TUJo{x) + YVoi^)<'^^ forallxGM^ 

Here, ivo{x) = —V x uo(x) = dyUo — dxVQ is the initial vorticity and r]o{x) := Ai — A2 is the (possibly 
complex-valued) spectral gap associated with the eigenvalues of gradient matrix Vuo(x). Moreover, 
these globally smooth solutions, ui(i, •), are 2iit -periodic in time. 

In [17] , Liu and Tadmor gave two different proofs of (j2.1b|) . One was based on the spectral dynamics 
of Aj(Vu); another, was based on the flow map associated with ()2.1ap . and here we note yet another 
version of the latter, based on the Riccati-type equation satisfied by the gradient matrix M =: Vui, 

Here {•}' := 5^ + ui • V denotes differentiation along the particle trajectories 
(2.2) Fo := {{x, t) I x{t) = ui{x{t),t), x{to) = xq}. 

Starting with Mq = M{to,xo), the solution of this equation along the corresponding trajectory Fq is 
given by 

M = e*-^/^ (/ + T-'^J (/ - e*-^/^) Mo) Mq, 
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and a straightforward calculation based on the Cayley-Hamilton Theorem (for computing the inverse 
of a matrix) shows that 



(2.3) niax |Vui| = rnax \M\ = rnax 



polynomial(r, e*"^/"^, Vuq) 



(1 - TWO - ^%) + 

Thus the critical threshold ()2.1bp follows. The periodicity of ui is proved upon integrating ui' = -Ju 
and x' = ui along particle trajectories Fq. It turns out both x{t) and ui(i, x(t)) are 2'kt periodic, 
which clearly implies that Ui(t,-) shares the same periodicity. It follows that there exists a critical 
Rossby number, Tc := rc(Vuo) such that the pressureless solution, ui(t, •), remains smooth for global 
time whenever r G (0,rc). This emphasizes the stabilization effect of the rotational forcing for a 
"large" class of sub-critical initial configurations, [17^ §1.2]. Observe that the critical threshold, Tc 
need not be small, and in fact, Tc = oo for rotational initial data such that r/g < 0, < \/— ??o- 
shall always limit ourselves, however, to a finite value of the critical threshold. Te- 
la the next corollary we show that in fact, the pressureless solution retains higher-order smoothness 
of the sub-critical initial data. To this end, we introduce the following notations. 

Notations. Here and below, || • ||m denotes the usual iif^'-Sobolev norm over the 2D torus and | • |oo 
denotes the L°° norm. We abbreviate a ^^y^ b for a < cb whenever the constant c only depends on the 
dimension m. We let Cq denote m-dependent constants that have possible nonlinear dependence on 
the initial data |/io|oo and |Vuo|oo- The constant, Co := Co • ||(/iO) uo)||mi will be used for estimates 
involving Sobolev regularity, emphasizing that Co depends linearly on the H'^-size of initial data, Hq 
and uo)||, and possibly nonlinearly on their L°°-size. 

Corollary 2.1. Fix an integer m > 2 and consider the pressureless system \2.1a\) subject to sub- 
critical initial data, uq £ H"^. Then, there exists a critical value Tc := Tc (Vuq) < oo such that for 
T £ (0, Tc] we have, uniformly in time, 

(2.4a) |Vui(t,-)|oo < Co, 

(2.4b) l|ui(t,-)llm < Co. 



Proof. We recall the expression for |Vui|oo in (|2.3p . By continuity argument, there exists a value t^ 
such that 1 — ruo ~ "y^o ^ ^/^ fo^' ah r G (0,rc) , which in turn implies (|2.4ap with a constant Co 
that depends on |Vuo|oo and Tc which also replies on the pointwise value of Vuo. 

Having control on the L°° norm of Vui, we employ the standard energy method to obtain the 
inequality, 

^l|ui(*r)l|m |Vui(t,-)|L-=||ui(t, Ollm- 

Since ui(t, •) is 27rr-periodic, it suffices to consider its energy growth over < t < 2irT < 2'kTc. 
Combining with estimate (j2.4ap and solving the above Gronwall inequality, we prove the estimate 
(l2lbl) . □ 

3. Second approximation - the linearized system 

Once we established the global properties of the pressureless velocity ui, it can be used as the 
starting point for second iteration of (jl.4p . We begin with the approximate height, /12, governed by 
(fLiip . 

(3.1) 5t/i2 + ui • V/12 + Q + /i2] V • ui = 0, h2{Q,-) = ho{-). 
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Recall that ui is the solution of the pressureless system ()2.1ap subject to sub-critical initial data uq, 
so that ui(t, •) is smooth, 27rr-periodic in time. The following key lemma shows that the periodicity 
of Ui imposes the same periodicity on passive scalars transported by such Ui's. 

Lemma 3.1. Let scalar function w be governed by 

(3.2) dtw + V ■ (uiw) = 

where ui{t,-) is a globally smooth, 2'itt -periodic solution of the pressureless equations Ii2.1a\) . Then 
w{t, •) is also 2'7TT -periodic. 

Proof. Let (j) := V x ui + denote the so-called relative vorticity. By (j2.1ap it satisfies the same 
equation w does, namely, 

dt^ + v- (ui0) = 0. 

Coupled with ()3.2p . it is easy to verify that the ratio w/(p satisfies a transport equation 

which in turn implies that w/(j) remains constant along the trajectories Tq in (j2.2p . But (I2.1al) tells 
us that ui' = 7U1, yielding ui(t, = er-^uo{xo). We integrate to find, x{2t:t) = x(0), namely, the 
trajectories come back to their initial positions at t = 2'kt. Therefore 

w w 

--{2ttt,xq) = — (0, xo) for all xq's. 

(j) 

Since the above argument is time invariant, it implies that w/(j){t,-) is 27rr-periodic. The conclusion 
follows from the fact that ui(t, •) and thus </>(f, •) are 27rr-periodic. □ 

Equipped with this lemma we conclude the following. 
Theorem 3.1. Consider the mass equation i3. 1]) on a 2D torus, T^, linearized around the pressureless 



velocity field ui and subject to sub-critical initial data (/iojUq) G i/™(T^) with m > 5. It admits a 
globally smooth solutions, h2(t,-) £ H^^^{T'^) which is 2TTT-periodic in time, and the following upper 
bounds hold uniform in time, 

(3.3a) \h2{t,-)\oo<Co(l + - 

V a 

(3.3b) ||/i2(t,-)IU-i <Co(l + - 

V a 

Proof. Apply lemma [3TT] with w := + /12 to (j3.ip to conclude that /12 is also 27rr-periodic. We 
turn to the examine the regularity of /i2. First, its L°° bound (j3.3ap is studied using the L°° estimate 
for scalar transport equations which yields an inequality for \h2\00 = |^2(i, Oloo 

^|/i2|oo < |V • Ui|oo(o-"^ + |/l2|oo)- 

dt 

Combined with the L°° estimate of Vui in (j2.4ap . this Gronwall inequality implies 

|/i2|oo <e^°*|/io|oo + -(e^°*-l 
a V 

As before, due to the 27rT-periodicity of /i2 and the subcritical condition t < Tc, we can replace the 
first t on the right with Tc, the second t with 27rT, and (j3.3ap follows. 
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For the F"""^ estimate prsbjl . we use the energy method and the Gaghardo-Nirenberg inequahty 
to obtain a similar inequahty for |/i2|m— i = m— 1 ) 

d ,,, ,, 1^ I nr. II /^l |, I A II II 

^||^2||m-l I VUi|oo||n'2||m-l + I ~ + l"'2|oo I ||Ul||m- 

Applying the estimate on ui in (j2.4p and the estimate on /12 in (j3.3ap . we find the above inequality 
shares a similar form as the previous one. Thus the estimate ()3.3bp follows by the same periodicity and 
sub-criticality argument as for ()3.3ap . We note by passing the Zmear dependence of Cq on ||(/io,uo)||m- 

□ 

To continue with the second approximation, we turn to the approximate momentum equation ()1.4bp 
with j = 2, 

1 1 

(3.4) dtU2 + U2 • VU2 + -V/l2 - -U2 = 0. 

a T 

The following splitting approach will lead to a simplified linearization of (j3.4p which is "close" to 
p.4p and still maintains the nature of our methodology. The idea is to treat the nonlinear term and 
the pressure term in (|3.4p separately, resulting in two systems for v ~ U2 and v ~ U2, 

(3.5a) 9iV + vV-v-ijv = 0, 

r 

(3.5b) dt^ + -Vh2 - - Jv = 0, 

(T r 

subject to the same initial data v(0, •) = v(0, •) = uo(-). 

The first system (|3.5ap , ignoring the pressure term, is identified as the pressureless system (|2.ip and 
therefore is solved as 

V = Ui, 

while the second system ()3.5bp . ignoring the nonlinear advection term, is solved using the Duhamel's 
principle, 

v(t,.) = e*^/- (uo(t,-)- r^^V/i2(s,-)ds) 







JJ/t 



ft g~sj/cr 

uo(t,-) - / Vh2{t,-)ds 

Jo o- 



= e"/^uo{t, •) + - J(I - e'-^/^)Vh2it, ■). 
a 

Here, we make an approximation by replacing h2{s, •) with h2{t, ■) in the integrand, which introduces 
an error of order r, taking into account the 27rT period of h{t, ■). 

Now, synthesizing the two solutions listed above, we make a correction to v by replacing e^'^/'^ug with 
ui. This gives the very form of our approximate velocity field U2 (with tolerable abuse of notations) 

(3.7a) U2 := ui + - J(I - e*-^/")V/i2(t, •)• 

a 

A straightforward computation shows that this velocity field, U2, satisfies the following approximate 
momentum equation, 

(3.7b) dtU2 + ui • Vu2 + -Vh2 - -ua^ = R 

a T 



LONG TIME EXISTENCE OF THE RAPIDLY ROTATING EULER EQUATIONS 

where 



(3.7c) 



R := -J{I -e"/^){dt + ui-V)Vh2{t,-) 



a 



T 



(bydMD) = J{I-e 



a 



(Vui)TV/l2+V((- + /l2)V-U 



1 



a 



Combining Theorem 13.11 on /i2(i, •) with Gaghardo-Nirenberg inequahty, we arrive at the following 
corollary on periodicity and regularity of U2. 

Corollary 3.1. Consider the velocity field U2 in 1^. 7| j subject to sub-critical initial data (/io,uo) G 
H"'{T'^) with m > 5. Then, U2(t, •) is a Ittt -periodic in time, and the following upper bound, uniformly 
in time, holds, 

||U2 - Ui||m_2 < Co- ( 1 + - 

In particular, sificG W'^iWm — C'o foT subcriticdl wc concludG that U2(^, •) has the Sobolcv TCQiildvity^ 

( T 

U2 U-2 < Co IH V ^ 

\ a a'^ 

We close this section by noting that the second iteration led to an approximate RSW system 
linearized around the pressuerless velocity field, ui, (j3.ip . (|3.7p . which governs our improved, 27rr- 
periodic approximation, (/i2(t, U2(i, •)) e i^^^H^^) x ^^"'"^(T^). 

4. Long time existence of approximate periodic solutions 



4.1. The shallow-water equations. How close is {h2{t, •), U2(t, •)) to the exact solution {h{t, •), u{t, •))? 

Below we shall show that their distance, measured in if™^^(T^), does not exceed tt-;—- Thus, 

l-gCot^ 

for sufficiently small 5, the RSW solution (/i, u) is "approximately periodic" which in turn implies its 
long time stability. This is the content of our main result. 

Theorem 4.1. Consider the rotational shallow water (RSW) equations on a fixed 2D torus, 
(4.1a) + u • V/i + + /i^ V • u = 

(4.1b) Stu + u- Vu + iv/i- iju = 

a T 



? 



5 



subject to sub-critical initial data (/io,uo) G H^-iT ) with m > 5 and ao := min(l +(t/io(-)) > 0. Let 

T 
T2 



denote the ratio between the Rossby number t and the squared Froude number a, with subcritical 
T < ''"c(Vuo) so that 112. lb\} holds. Assume a < 1 for substantial amount of pressure forcing in (4.1b). 
Then, there exists a constant Cq, depending only on m, t^, ao and in particular depending linearly 
on ||(/io, uo)||m! such that the RSW equations admit a smooth, "approximate periodic" solution in the 
sense that there exists a near-by 2ttt -periodic solution, (/i2(t, •), U2(t, •)), such that 

(4.2) ||p(t,-)-P2(t,-)IU-3 + ||u(t,-)-U2(t,-)IU-3 < ^_^Cots - 

Here p is the "normalized height" such that 1 + ^ap = \/l + ah, and correspondingly, p2 satisfies 
1 + ^crp2 = VT+'ah^. 

R follows that the life span of the RSW solution, t < ts '■= ln{5~^) is prolonged due to the rapid 
rotation 5 <^ 1, and in particular, it tends to infinity when (5 — > 0. 
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Proof. We compare the solution of the RSW system (j4.1ap . (|4.1bp with the solution, (/i2;U2), of ap- 
proximate RSW system (]3.ip . (j3.7p . To this end, we rewrite the latter in the equivalent form, 

(4.3a) dth2 + U2-\/h2+ (^^ + h2^V ■U2 = (ua - ui) • V/i2 + + /i2^ V • (u2 - ui) 

(4.3b) C?tU2 + U2 • VU2 + — V/l2 - -'/U2 = (U2 - Ui) • VU2 + -R. 

a T 

The approximate system differs from the exact one, (j4.1ap . ()4.1bp . in the residuals on the RHS of 
(|4.3ap . (l4.3bp . We will show that they have an amplitude of order 6. In particular, the comparison in the 
rotationally dominant regime, (5 <C 1 leads to a long-time existence of a smooth RSW solution which 
remains "nearby" the time-periodic solution, (/12, U2). To show that (/i2, U2) is indeed an approximate 
solution for the RSW equations, we proceed as follows. 

We first symmetrize the both systems so that we can employ the standard energy method for 
nonlinear hyperbolic systems. To this end. We set the new variable ("normalized height") p such that 
1 + ^ap = \/l + ah. Compressing notations with U := (p, u)^, we transform (j4.1ap . (|4.1bp into the 
symmetric hyperbolic quasilinear system 

(4.4) 9fU + 5(U,VU) +i^[U] = 0. 

Here B(F, VG) := Ai{F)Gx + A2{F)Gy where Ai, A2 are bounded linear functions with values being 
symmetric matrices, and K\F] is a skew-symmetric linear operator so that (i^[F], F) = 0. By standard 
energy arguments, e.g. [12] , [13] , [E] ) , the symmetric form of (14. 4p yields an exact RSW solution U, 
which stays smooth for finite time t ^ 1. The essence of our main theorem is that for small (5's, 
rotation prolongs the life span of classical solutions up to t ~ 0(ln5^^). To this end, we symmetrize 
the approximate system (|4.3ap . ()4.3bp . using a new variable p2 such that 1 + |crp2 = \/l + (t/i2- 
Compressing notation with U2 := (p2,U2)^, we have 

(4.5) 9tU2 + S(U2,VU2) + i^(U2) = R 

where the residual R is given by 

(U2 - ui) • Vp2 + (I + P2) V • (U2 - Ui) 
(U2 - Ui) • VU2 - R 



R 



with R defined in (j3.7cp . We will show R is small which in turn, using the symmetry of (j4.4p and 
(|4.5p . will imply that ||U — U2||m-3 is equally small. Indeed, thanks to the fact that H'^~^(T'^) is an 
algebra for m > 5, every term in the above expression is upper-bounded in H"^~'^, by the quadratic 
products of the terms, ||ui||m, ||P2||m-i5 ||u2||m-2, ||u2 — ui||m-2, up to a factor of 0(1 + ^). The 
Sobolev regularity of these terms, ui, U2 and p2 is guaranteed, respectively, in corollary 12.11 corollary 
13.11 and theorem 13.11 Moreover, the non- vacuum condition, 1 + (j/iq > ckq > 0, implies that 1 + ah2 
remains uniformly bounded from below, and by standard arguments (carried out in Appendix A), 
||P2||m-2 < Co{l + t/ct). Summing up, the residual R does not exceed. 



(4.6) ||R|U_3 <Cl{S + l + ... + ^] <Cl6, 



a a' 



for sub-critical r G (0, Tc) and under scaling assumptions 5 < 1, cj < 1. 

We now claim that the same C'((^)-upperbound holds for the error E := U2 — U, for a long time, 
t ^ ts. Indeed, subtracting (j4.4p from (j4.5p . we find the error equation 



dtE + B{E, VE) + K[E] = -5(U2, VE) - B{E, VU2) + R. 
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By the standard energy method using integration by parts and Sobolev inequaUties while utihzing the 
symmetric structure of B and the skew-symmetry of K, we arrive at 

^l|E||m-3 ~m l|E||^_3 + ||U2 ||m-2 ||E||^_3 + ||R||m_3 ||E||m_3. 

Using the regularity estimates of U2 = (p2,U2)^ and the upper bounds on R in (14. 6p . we end up 
with an energy inequality for ||E(t, •)||m-3, 

^||E||m-3 ~m l|E||m~3 + C'o||E||m--3 + Cq^, I|E(0, •)l|m-3 = 0. 

A straightforward integration of this forced Riccati equation (consult for example, [16, §5]), shows 
that the error ||E||m-3 does not exceed 

(4.7) ||U(t,-)-U2(t,-)IU-3 < T 



In particular, the RSW equations admits an "approximate periodic" ff™~'^(T^)-smooth solutions for 
t < ^ln{dde-^) iov 6 ^1. □ 

Remark 4.1. The estimate on the actual height function h follows by applying the Gagliardo-Nirenberg 
inequality to h - h2 = p{l + jp) -^2(1 + JP2) = (p-p2)(l + j{p - P2) + fP2), 

\\h{t,-)-h2{t,-)\\m-3^ ^^_^Cotsy 

Our result is closely related to observations of the so called "near-inertial oscillation" (NIO) in 
oceanography (e.g. [22]). These NIOs are mostly seen after a storm blows over the oceans. They 
exhibit almost periodic dynamics with a period consistent with the Coriolis force and stay stable for 
about 20 days which is a long time scale relative to many oceanic processes such as the storm itself. 
This observation agrees with our theoretical result regarding the stability and periodicity of RSW 

solutions. In terms of physical scales, our rotationally dominant condition, 6 = Yjjj ^ ^' Pi'^'^i^^s ^ 

physical characterization of this phenomenon. Indeed, NIOs are triggered when storms pass by (large 
U) and only a thin layer of the oceans is reactive (small aspect ratio H/L). Upon using the multi-layer 
model ([18j §6.16]), we consider scales / = 10~^s~^,L = 10^m,i? = lO^m, C/ = lms~^,g = O.Olms"^ 
(reduced gravity due to density stratification ~ consult §1-3]). With this parameter setting, 
5 = 0.1, and theorem 14.11 implies the existence of smooth, approximate periodic solution over time 
scale ln(5^^)L/C/ ~ 2 days. We note in passing that most cyclonic storms on the Northern Hemisphere 
rotates clockwise, yielding a negative vorticity, ojq = dyUQ — dxVo < 0, which is a preferred scenario of 
the sub-critical condition (|2.1bp assumed in theorem 14. 1[ 

4.2. The isentropic gasdynamics. In this section we extend theorem 14.11 to rotational 2D Euler 
equations for isentropic gas, 

(4.8a) dtp + V-{pu) = 

(4.8b) (9tu u • Vu p^^Vpip) - /u-L = 0. 

Here, u := {u^^\u^'^'^)~^ is the velocity field, p is the density and p = p{p) is the pressure which for 
simplicity, is taken to be that of a polytropic gas, given by the 7-power law, 

(4.8c) p{p) = Ap'^. 
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The particular case A = 5/2,7 = 2, corresponds to the RSW equations (|l.lap . (|l.lbp . The following 
argument for long term existence of the 2D rapidly rotating isentropic equations applies, with minor 
modifications, to the more general pressure laws, p{p), which induce the hyperbolicity of (j4.8ap . 
We first transform the isentropic Euler equations (j4.8ap into their nondimensional form, 

dtp + u-Vp+ (^^ + p^V -u = 

dtu + u-Vu + ^V(l+apy-^ --Ju = 
cr^ T 

where the Mach number a plays the same role as the Froude number in the RSW equation. In order 
to utilize the technique developed in the previous section, we introduce a new variable h by setting 
1 + ah = {1 + apY^^, so that the new variables, (/i, u), satisfy 

(4.10a) + u • V/i + (7 - 1) + /i^ V • u = 0, 

(4.10b) Sju + u- Vu + iv/i- -Ju = 0. 

a T 

This is an analog to the RSW equations (j4.1ap . ()4.1bp except for the additional factor (7 — 1) in the 
mass equation (|4.10ap . We can therefore duplicate the steps which led to theorem 14.11 to obtain a long 
time existence for the rotational Euler equations (j4.10ap . ()4.10bp . We proceed as follows. 

An approximate solution is constructed in two steps. First, we use the 27rr-periodic pressureless 
solution, {h\ = constant, ui(t, •)) for sub-critical initial data, (/io,uo). Second, we construct a 27rT- 
periodic solution (/i2(i, •)> U2(i')) ^ the solution to an approximate system of the isentropic equations, 
linearized around the pressureless velocity ui, 

dth2 + ui • V/12 + (7 - 1) (^^ + V • ui = 0, 

U2 := ui + ^ J - e*^/") V/i2(t, •)• 

In the final step, we compare (/i, u) with the 27rr-periodic approximate solution, (/12, U2). To this end, 
we symmetrize the corresponding systems using U = (p, u)"*^ with the normalized density function 

ip satisfying 1 + \\j t^^'^V = \/l + crh. Similarly, the approximate system is symmetrized with the 
variables U2 = (^2, U2) where 1 + \ 7:pi(^V2 = \/l + <^^2- We conclude 



Theorem 4.2. Consider the rotational isentropic equations on a fixed 2D torus, \4.9a% 1^4. 9a ), subject 

to sub-critical initial data {pQ, uq) G i?™(T^) with m > 5 and := min(l + crpo(")) > 0- 

Let 

denote the ratio between the Rossby and the squared Mach numbers, with sub- critical r < Tc (Vuq) 
so that 112. lb]) holds. Assume a < 1 for substantial amount of pressure in ^.9a ). Then, there exists 



a constant Co, depending only on m, ||(poi uo)||m, Tc o,nd ao, such that the RSW equations admit a 
smooth, "approximate periodic" solution in the sense that there exists a near-by 2ttt -periodic solution, 
{p2{t, ■),U2{t, •)) such that 

(4.11) Mt, •) -P2it, ■)\\m-3 + ||u(i, •) - U2(t, OIk-S < Y^^' 
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7 — 1 

Here, p is the normalized density function satisfying 1 + ap = (1 + (jp)^~ , and p2 results from the 
same normalization for p2- 

It follows that the life span of the isentropic solution, t ^ ts ■= 1 + ln((5~^) is prolonged due to the 
rapid rotation 5 <C 1, and in particular, it tends to infinity when 5^0. 

Remark 4.2. For the actual density functions, p — P2 = ^[(1 + crp)''-'^ — (1 + 0-^2)^"^] = Jq C^[1 + 

(Tie{p-P2)+P2)]^'^de 

gCot 

\\p{t, •) - P2it, OIU-S ^ 

in the physically relevant regime 7 G (1,3). 

4.3. The ideal gasdynamics. We turn our attention to the full Euler equations in the 2D torus, 

dtp + v-{pu) = 0, 

dtu + u-Vu + p-^Vp{p,S) = fJu, 
dtS + u-VS = 0, 

where the pressure law is given as a function of the density, p and the specific entropy S, p{p, S) := 
p'^e^ . It can be symmetrized by defining a new variable - the "normalized" pressure function, 

p := -^p 2t , 
7-1 

and by replacing the density equation ()4.12ap by a (normalized) pressure equation, so that the above 
system is recast into an equivalent and symmetric form, e.g., |12j.[llj 

e^dtp + e'^u • Vp + C^e^pV • u = 0, 

dtu + u-Vu + C^e^pVp = fJu, := 

dtS + u-VS = 0. 

It is the exponential function, e'^, involved in triple products such as e^pVp, that makes the ideal gas 
system a nontrivial generalization of the RSW and isentropic gas equations. 
We then proceed to the nondimensional form by substitution, 

u^Um', p^P{l + ap'), S = ln{pp-^) ^ln{PR-^) +aS' 

After discarding all the primes, we arrive at a nondimensional system 

(4.13a) e'^^dtp + e'^^u -Vp + C^ (^^-^ + e^^p) V • u = -C^^V ■ u, 

(4.13b) dtu + u-Vu + Cy( ^ ^ ~ + e'^^p] Vp = -C^-Vp+-Ju, 

\ a J a T 

(4.13c) dtS + u-VS = 0, 

where a and r are respectively, the Mach and the Rossby numbers. With abbreviated notation, 
U := {p, u, S)~^ , the equations above amount to a symmetric hyperbolic system written in the compact 
form, 

(4.14) AoiS)dtU + Ai(U)9,U + A2{V)dyV = K[V]. 
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Here, Ai{i = 0,1,2) are symmetric-matrix-valued functions, nonlinear in U and in particular Aq is 
always positive definite. The linear operator K is skew-symmetric so that (i^[U],U) = 0. 

Two successive approximations are then constructed based on the iterations (jl.4p . starting with 
J = l, 

Pi = constant, 
9tui + ui • Vui = -Jui, 

T 

Si = constant. 

Identified as the pressureless solution, ui is used to linearize the system, resulting in the following 
approximation 

(4.15a) 9tP2 + ui • Vp2 + C^P2V • U2 = -C^^V • ua, 

(4.15b) U2-U1 = Ij(/-e*-^/nC7^e"^2(l + ap2)Vp2, 

(4.15c) dtS2 + ui-VS2 = 

The 27rr-periodicity and global regularity of U2 := (pa, U2, 52)^ follow along the same lines outlined 
for the RSW equations in section [3] (and therefore omitted), together with the following nonlinear 
estimate for e'^^ , 



g {asy 



2. n w^Wm = ^ I — O-'S 



for the latter, we apply recursively the Gagliardo-Nirenberg inequality to typical terms ||(cj5)-' 
Notice the entropy variable (both the exact and approximate ones) always satisfies a transport equation 
and therefore is conserved along particle trajectories, which implies that the norm of the entropy 
variable is an invariant. Thus, we arrive at an estimate 

(4.16) lle-^^-lIU <(Tao||5|U. 

Of course, the same type of estimate holds for the approximate entropy, ^2. 

Finally, we subtract the approximate system (j4.15p from the exact system ()4.14p . arriving at an 
error equation for E := U — U2 that shares the form as for the RSW system in Section [4.H except 
that Ai{\J) — Ai{\J2) / Ai{\J — U2) due to nonlinearity which is essentially quadratic in the sense 
that, 

Pi(U)-Ai(U2)|U < ||U-U2||^ + ||U-U2|U, i = 0,l,2, 

Pi(U) - Ai(U2)||w'i,oo < ||U-U2||^i.oo + ||U-U2||vKi.-, i = 0,l,2. 
where n > 2. This additional nonlinearity manifests itself as three more multiplications in the energy 
inequality, 

^||E|U_3 < ||E||^_3 + ... + ||E|U_3 + 5, ||E(0, Ollm-a = 0, 

whose solution (- developed around a simple root of the quintic polynomial on the right), has the 
same asymptotic behavior as for the quadratic Riccati equations derived in the previous sections. 



^Consider a typical term of Ai, e.g. e'^^p. Applying (|4.16|l together with Gagliardo-Nirenberg inequality to e'^'^ — 
g<TS2 ^ g,TS2(g<T(s-S2)-i)^ ^Yiow We"'^ - e-'^^lln < \\S - S2||n. The estimate on ||e'"^p - e'"'^2p2 ||n then follows by 

applying identity 06 — 02&2 = (o — a2){b — 62) + (a — 02)^2 + a2(6 — ^2) together with the triangle inequality and the G-N 
inequality. Here regularity of S2 and P2 is a priori known. 
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Theorem 4.3. Consider the (symmetrized) rotational Euler equations on a fixed 2D torus 

subject to sub-critical initial data {pq^vlq, Sq) G //'"(T^) with m> 5. 

Let 

denote the ratio between the Rossby and the squared Mach numbers, with subcritical r < Tc{nablauQ) 
so that 112. lb\} holds. Assume a < 1 for substantial amount of pressure forcing in (4.13b\ ). Then, 
there exists a constant C, depending only on m, \\{po,uo, So)\\m, Tc such that the ideal gas equations 
admit a smooth, "approximate periodic" solution in the sense that there exists a near-by 2'KT-periodic 
solution, {p2{t, ■),U2{t, ■), S2{t, ■)) such that 

(4.17) Mt, •) -p2it, OIU-s + ||u(t, •) - Mt, ■)\\m-3 + \\Sit, •) - S2{t, OlU-a < ^_^Cots - 

It follows that the life span of the ideal gas solution, t ts ■= ln{5^^) is prolonged due to the rapid 
rotation 5^1, and in particular, it tends to infinity when 5^0. 

5. Appendix. Staying away from vacuum 
We will show the following proposition on the new variable p2 defined in section 14.11 
Proposition 5.1. Let p2 satisfies 

(5.1) 1 + ^ap2 = \/l + ah2 



where /i2 is defined as in liS.l]) . that is, 

(5.2) at/i2 + ui • V/l2 + (^^ + /i2^ V • ui = 

subject to initial data /i2(0, •) = /io(") that satisfies the non-vacuum condition 1 + (t/io(") > ckq > 0. 
Then, 

\P2\00 < Co (1 + 

\\p2\\n < Co (1 + - 

\ a 

The proof of this proposition follows two steps. First, we show that the L°° and norms of 
^2(0,-) dominated by /i2(0, •) due to the non-vacuum condition. Second, we derive the equation 
for p2 and obtain regularity estimates using similar techniques from section 14.11 

Step 1. For simplicity, we use p := ^2(0, •) and h := /i2(0, •). 

Solving ()5.1|) and differentiation yield 

2h Vh 
P = ■■ ) V p 



Vl + crh + l' yjl + ah' 

Clearly, \p\oo ^ |^|oo- The above identities, together with the non- vacuum condition imply 

IIpIIi < 2\\hh and \Vp\l^ < 

For higher derivatives of p, we use the following recursive relation. Rewrite ()5.ip as, p + \ap^ = h 
and then take the A:-th derivative on both sides 



D^p + ^a2pD^p + ^cj (p^{q^) - 2pD^p^ = D' 



h 
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SO that taking L norm of this equation yields 



(1 + ]^ap)D^p 



1 



D''{q^) - 2pD''p < \\D''h\\o 







Furthermore, we find / > y^ao||-D'^p||o by (|5.ip and the non-vacuum condition. We also find II 
|Vp|oo||p|||fc|-i by Gagliardo-Nirenberg inequalities. Thus we arrive at a recursive relation 

||p|||fe| <Co(||p|||fc|-i + ||/i|||fc|) 

which implies that the norm of P2(0, •) — p is dominated by ||/i2(0) Oil™ ~ 
Step 2. We derive an equation for p2 using relation (jS.ip and equation (j5.2p . 



dtP2 + 2ui • Vp2 + + V • ui = 0. 

This equation resembles the formality of the approximate mass equation (j3.ip for /i2 and thus we 
apply similar technique to arrive at the same regularity estimate for p2, 

|p2(t,-)|oo < Co(l + - 

\ a 
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